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Abstract 

The cluster state model for quantum computation [Phys. Rev. Lett. 
86, 5188] outlines a scheme that allows one to use measurement on a large 
set of entangled quantum systems in what is known as a cluster state to 
undertake quantum computations. The model itself and many works ded- 
icated to it involve using entangled qubits. In this paper we consider the 
issue of using entangled qudits instead. We present a complete frame- 
work for cluster state quantum computation using qudits, which not only 
contains the features of the original qubit model but also contains the 
new idea of adaptive computation: via a change in the classical compu- 
tation that helps to correct the errors that are inherent in the model, the 
implemented quantum computation can be changed. This feature arises 
through the extra degrees of freedom that appear when using qudits. Fi- 
nally, for prime dimensions, we give a very explicit description of the 
model, making use of mutually unbiased bases. 



1 Introduction 

The problem of building a quantum computer that has computational power 
greater than that of classical computers is one that has concerned both theo- 
rists and experimentalists for a number of years. The recent development of 
the cluster state model for quantum computation El El j m which remark- 
ably quantum computation is achieved essentially by tailored measurement of 
entangled states of a large number of qubits, has not only lead to significant 
implications in our understanding of quantum computation and quantum dy- 
namics, but also to new experimental approaches to the possible realisation of 
large scale quantum computers pfl 1221 1?)]. 

Although much study has been devoted to this model El El El EH] 
this work has almost exclusively been concerned with using qubits as the bàsic 
physical resource. However, many quantum systems cannot be treated as simple 
two-level systems, but rather as multi-level systems. Furthermore, recent work 



(i has suggested that qutrit based quantum computation schemes yield in some 
sense the most emcient implementation of quantum computation. 

In this paper we are going to present a generalisation of the cluster state 
model, by replacing the cluster state of qubits with more general cluster states 
of qudits. Such a generalisation has already been presented in but in this 
paper we will take a much more explicit approach, in which we try and present 
the cluster state model by presenting each constituent of the model in as simple 
a way possible. This allows us to see the critical ingredients that make cluster 
state QC possible. 

The paper will be laid out as follows. In section |21we will present a gen- 
eralisaton of an idea used in [2] which we will call one dit teleportation. This 
forms the core to our approach to cluster state QC. One dit teleportation takes 
the form of a simple two-qudit circuit identity which can be thought of as im- 
plementing a certain quantum gate (namely the d-dimensional analogue of the 
Fourier transform gate) by measurement. We then proceed to show that this 
identity can be used to implcment a much wider class of gates, and how this pro- 
cess can be turned into a quantum computational paradigm, which essentially 
involves a series of measurements on a entangled state of qudits to implement 
quantum computation. While this approach turns out to not be the best ap- 
proach for some types of cluster states, it ties in our cluster state model for QC 
heavily with the circuit model, henec giving us an alternative and (in the eyes 
of the author) a simpler way of looking at the model. 

In section|31we will first discuss the issue of adapting measurements based on 
previous measurement results to allow us to implement any quantum operation 
we wish. Measurements on a quantum system by their very nature are statis- 
tical; and in the cluster state paradigm this randomness can be thought of as 
introducing errors into our computation which can be corrected for using clas- 
sical computation. This notion is known as the issue of adaptive measurement 
and is crucial in the theory of cluster state QC. We will also introduce a new 
phenomena which we call adaptive computation, a feature that will be unique 
to cluster state QC on qudits, arising due to the extra degrees of freedom we 
are presented with. The bàsic idea is that by changing the method by which we 
correct the errors associated with a measurement outeome, we can implement 
different quantum gates. We will discuss this in further detail at the appropriatc 
time. We conclude this section with an example of the simplicity of the one-dit 
teleportation approach which also makes use of adaptive computation. 

In section0J we briefly review the stabiliser method for cluster states, which 
is the most common approach in most of the existing literature. For details 
of the model we refer the reader to the appropriate paper, but we do give a 
comparison of the two approaches, highlighting the fact that both approaches 
are indeed useful. 

In section we will discuss the problem of generating a universal set of 
quantum gates for quantum computation, giving a general method by which 
quantum gates can be realised, and then discussing a connection of this method 
with the theory of mutually unbiased bases. While this connection is not a 
fundamental one, using the known theory for MUBs gives us a very elegant 
method of establishing measurement patterns for implementing a universal set 
of quantum gates. We will use the connection to give explicit calculations and 
details of the measurements required to generate certain quantum gates in prime 
dimensions. We conclude with a discussion in the final section. 
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2 One dit teleportation and cluster state quan- 
tum computation 



2.1 Preliminaries: a bàsic measurement result 

For clarity, we will start by giving a very bàsic result. We assume throughout this 
paper that all Hilbert space(s) involved here are of dimension d. Let 
be a Standard (computational) basis for our Hilbert space. We define measuring 
in the basis defined by U for a unitary U £ U (d) to mean measurement in the 
orthonormal basis {U\k)}^Z^ (i.e. the vectors defined by the columns of U 
when U is described in matrix form in the computational basis). Our result is 
as follows: 

Lemma 1. Given a bipartite state \^}, measuring one of the systems in the 
basis defined by U and discarding is equivalent to applying to the first system, 
measuring it in the computational basis and then discarding it. 

Proof. Let = J2k^\^) <8> \ipk)- Measuring the first system in the basis 
defined by U leaves the second system in the state \ipk) / \J (V'felV'fc) with prob- 
ability (ipkl^k)- Applying W to the first system and then measuring it in the 
computational basis can easily be seen to have the same effect. 

This lemma gives us an alternative way to generate the appropriate statistics 
for measurement of one system of a bipartite system in a given basis. 

2.2 One dit teleportation 

The idea of one dit teleportation is a many-level generalisation of one-bit tele- 
portation that was first presented in We will first simply present this 
relatively simple idea, and then show how it is integral to cluster state quantum 
compuation. 

Wc first need some definitions. The Fourier transform basis is defined by 

d-1 

!+,-> = XX* l*> (i) 

fc=0 

where u = e 27rí / d , the primitive dth root of unity. We will often denote |+o) 
by simply |+). We also need to introduce a set of operators known as the 
generalised Pauli operators |14j : 

Z = |> fe |fc)(fc|, X = J2\k-l)(k\ (2) 

k=0 k=0 

where we are using modulo d arithmetic within the bras and kets. Finally, we 
define the Fourier gate F by 

d-1 d-1 

F = -r d H ^|j>(fc| = £|+*><fc|. (3) 

v j.k—O k—O 

This gate is the d-dimensional Quantum Fourier Transform gate, and in d = 2 
is equivalent to the Hadamard gate. 
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Figure 1: The circuit diagram for one-dit teleportation. 

The notion of one-dit teleportation is given by the circuit identity in Figure 
^ where we are using the d dimensional analog of controlled gates: 

d-1 d-1 

CZ = Y, \k){k\ ® Z k = e 27Ttkl/d \k)(k\ ® \l)(l\ (4) 

fc=0 k,l=0 

and the mctcr rcpresents measurement in the computational basis, with outcome 
k corresponding to the state \k). Furthermore, each of the possible outcomes 
occurs with equal probability 1/d. 

Proof. First we note two important facts, both of which can be deduced by 
elementary means: 

Z 3 l+fc) = |+j+fc) (modulo d addition in ket index); (5) 

X*\+ k )=wt k \+ k ). (6) 
Let = J2k a k\k)- Then 

1^)1+) = 5> fc |*>|+) 

fc 

J2a k \k)\+ k ) 

k 

E«fc l + fc) l + fc) ^1*) 

fc 

where we calculate the effect of the CZ gate using equation Now, we can 
write l'S) as 

= -^E^E a ^ fc i+fc) 

V 3 k 

= 4E^E a ^^i fc > 

3 fc 

= -^Ei^ F ^> 

3 

where the penultimate line follows from equation JHJ above; hence, when mea- 
suring the first qudit, an outcome j (occuring with probability yields a 
state X 3 F \ip) for the second qudit. 



fc 

cz 



k 
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2.3 Application to cluster state quantum computation 



One way we can think about the above circuit identity is that it gives us a very 
simple scheme for implementing the quantum Fourier gate F: Given a qudit 
pair in states \ip) , |+) entangled with a controlled-Z interaction, if we measure 
the first (1^)) qudit in the basis, and obtain the measurement outcome j, the 
second qudit ends up in the state X 3 F \tp). The X 3 operator can be thought of 
as an error operator, and we will discuss this shortly. 

Furthermore, from equation we note that the controlled-Z gate is sym- 
metric in the two qudits, and so any phase transformation in the computational 
basis (that is, a quantum gate of the form 

Z(a)=X)e < °*|A)(*l- 
fe 

with a £ [0, 2ir] d ) that acts on the first (upper) qudit commutes through the 
controlled-Z gate. This means that measuring the first system in the basis 
defined by {FZ(a))\ which (by lemma^l is the same as applying FZ(&) to 
the system before measurement in the computational basis, is equivalent to 
teleporting the initial state Z(&) \ïp). This leads to the more general circuit 
identity given in Figure[2l meaning we can implement the quantum gates FZ(&) 
using this method. 



l+> -z 



Z{a) - F - 



X m FZ(&) \ijj) 



Figure 2: The more general circuit diagram for one-dit teleportation. 

If the transformation \ip) — > FZ(a) \tp) forms our entire quantum compu- 
tation, then since the powers of X merely permute the computational basis 
elements, and all known quantum algorithms conclude with a measurement in 
this basis, we can simply correct for the Pauli error X 3 via a classical compua- 
tion. 

This is the simplest possible example of cluster state computation: We take 
two qudits, one of which is our initial state and an output qudit, which 
is initialised in the state |+). After entangling the pair using a controllcd-Z 
interaction, and an appropriate measurement, we can map our input state 
onto an output state FZ(a) \ip), up to a Pauli error X 3 . 

If we wished to apply further quantum operations to our state, we could do 
so by using a linear cluster of qudits entangled in this way (Figure |3J, and then 
measuring along the cluster appropriately to implement the appropriate prod- 
uct of gates. Since measurement on a qudit and an interaction between another 
pair of qudits commute, this is equivalent to several one-dit teleportations, with 
the output from one teleportation becoming the input to the next (see Figure 
|2J|. The measurements after the first however may have to be adapted due to 
the extra X factors brought in by previous measurements, so that the desired 
quantum operation is implemented. This issue of compensating for the random- 
ness of previous measurement outeomes by a change of subsequent measurement 
bases is known as adaptive measurement, and will be discussed in section 3. For 
now however, we will assume that this issue is surmountable. 
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These ideas lead us to our dcfinition of the cluster states we are going to use 
for computation: 

1. We prepare a set of qudits, each in the state |+) (except perhaps for an 
'initial' set of qudits which are in an potentially entangled state |^>); this 
represents the input to the algorithm); 1 

2. Finish creating the cluster by interacting 'neighbouring' pairs of qudits 
via a controlled-Z gate. 

For qubits, this is identical to the original framework for cluster state QC P3G]- 
Cluster state computation then takes the form of a series of measurements 
on specified qudits, along with a scheme on how to adapt the measurements 
based on the previous measurement results. This definition also gives a possi- 
ble mcthod for generating these cluster states: suppose we have d-dimensional 
physical systems that we can prepare in the state and we can allow them to 
interact with an appropriate interaction Hamiltonian between the systems e.g. 

H I = -gJ2kl\k)(k\ ® | Z> </ 1 

k,l 

where g is some coupling constant. If we let this interaction run for a time 
t = 2tt/ gdh, then the time evolution operator for the two qudits takes the form 
of a controlled-Z gate. Physical systems which could be used to implement 
this scheme are discussed in JH] and appropriate references therein. For the 
remainder of this paper, we are going to discuss the mathematical framework 
of this model. 

Before we continue, we make one final remark about our cluster states: they 
satisfy a property known as maximal connectedness, which for our purposes we 
can state as follows: Given an arbitrary cluster state, if we measure any qudit 
of the cluster in the computational (Z) basis, the remaining states are left in 
a cluster state, up to some Pauli Z errors. For qubit clusters, this issue is 
first discussed in \T2\ . This fact can be used to design specific clusters from 
larger base clusters (e.g. a grid) as by the above, Z-measurements on a qudit 
effectively remove it from the cluster state. This is discussed in more detail in 
appendix E] 

Suppose we can use linear clusters to implement any single qudit operation. 
We can then use more general multi-dimensional clusters to allow the implc- 
mentation of more complicated multi-qudit algorithms. One possible model is 
where rows of linear clusters, each representing a logical qudit (i.e. correspond- 
ing to a single 'qudit wire' in the circuit model) are used to implement single 

1 Alternatively, some initial measurements on a cluster state where before any interactions 
are applied all qudits are initially in the state |+) could be used on a larger cluster to prepare 
this initial state. 



Figure 3: A linear cluster of 6 qudits. Each dot represents a physical qudit, and 
the lines an interaction between the two adjoined qudits. 
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Figure 4: A quantum circuit for two successive one-dit teleportations. The 
circuit diagram is arranged to illustrate that the teleportation between the upper 
two qudits is independent of any further structure in the cluster (in this case a 
third qudit entangled to the second). The final state will take the form \ip') — 
X a(m 1 ,m 2 ) Z 6(m 1 ,m 2 ) FZ ( c ( b; mi ^FZ(a) \ip) , where c is a function oïm 1 because 
of the generalised Pauli error introduced by the first measurement may require 
the second measurement to be adapted. 



qudit operations, and entanglcment (again in the form of a controlled-Z inter- 
action) between the rows acts as an interaction between logical qudits. We will 
use this topology for cluster states a number of times throughout this paper; 
however, there are also a number of other ways of establishing this interaction 
between logical qudits. 2 It is a known fact that the set of one qudit gates along 
with any interacting two-qudit gate is sufficient to create any unitary gate on 
an arbitrary number of qudits (i.e. this set of gates is a universal set for multi- 
qudit quantum gates) [Uj- Hence by establishing single-qudit universality, we 
obtain universality for an arbitrary number of qubits on an appropriate cluster. 
A diagram to illustrate this point is given in Figure 



Figure 5: An example of a topology for cluster states that can be used to 
establish universal quantum computation. The rows represent a logical qudit, 
with one-dit teleportation the mechanism for implementing single-qudit gates. 
The interactions between rows are used to implement two-qudit gates, and so 
we can implement any arbitrary multi-qudit unitary operation. 

We note that although in this paper we study single-qudit universality on 
linear clusters, since ultimately these ideas are to be applied to larger multi- 
dimensional clusters, we are not contradicting the theorem of Nielsen |16j that 
says that quantum computation on linear clusters can be efficiently simulated on 
a classical computer (although the proof is for qubits, the same proof applies to 
qudits), as we ultimately want to use multi-dimensional clusters for multi-qudit 
QC. 

For qubits, the set of gates given by C2 = {FZ(&) | a E [0,27r] 2 } (i.e. those 
that can be implemented via one-&zí teleportation) is a universal set for single- 
qudit quantum gates. 3 In sectionElwe will consider the same question for the 
set C d = {FZ(a) | ae [0,27r] rf }. 

2 The computational power of cluster states of many shapes are studied in |13| . 

3 This is equivalent to the well known Euler decomposition for 2-dimensional unitàries. 
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Finally, we note that our defmition of qudit cluster states is the same as that 
used in [10] . In that paper the stabiliser formalism for cluster state computation 
(introduced in and discussed later in section 0} is heavily used. We will use 
a much more direct approach to establish the theory. 

3 Adaptive measurements and adaptive compu- 
tation 

In this section we outline how the unwanted X m (generalised Pauli) operators 
that appear in our cluster state model can be corrected using classical computa- 
tion, and by adapting, based on previous measurement results, the bases furthcr 
measurements are taken in. This idea is known as adaptive measurements. We 
will also show that by changing our correction method for these errors, we can 
extend one-dit teleportation to a further class of quantum gates; we call this 
adaptive computation. 

3.1 Adaptive measurements for single qudit operations 

In using linear clusters and successive one-dit teleportation defined in the pre- 
vious section to implement one-qudit quantum gates, generalised Pauli errors of 
the form X m are introduced to the overall state after each measurement. What 
we aim to show here is a remarkable feature of the original cluster state model 
that spreads to our generalisation, that we can compensate for these errors using 
solely classical computation. 

Suppose we have a state \ip) ~ X X Z Z U \<f>). The U represents the quantum 
gate applied so far, and the powers of X and Z are the generalised Pauli errors 
that we have due to previous measurements. Suppose now we apply one-dit 
teleportation to the pair \ip) \+) entangled by a controlled-Z operation as usual, 
and we obtain the new state X m FZ(a) The final state can be written as 
\ip') = X rn FZ(a)X x Z z U |0) .The problem that we appear to have is that we are 
implementing the gate FZ(a)X x Z Z U rather than the intended FZ(a)U. What 
we want is to move the powers of X and Z through the F and Z(a) so that 
they appear on the left of the overall gate FZ(&)U, so we can treat the Pauli 
operator as an error operator. To do this, we can make use of the following 
identities: 

Z(a)X = XZ(a'); « = a fe _ x ) (7) 

Z(a)Z = ZZ(a) (8) 
FZ = XF : FX = Z~ X F (9) 

and hence (ignoring any changes in the irrelevant overall phase), 

= X m FX x Z z Z(a íx) )U\<f>) 
= X m Z- x X z FZ(aW)U\<f>) 

= uj xz X m+z Z' x FZ{a {x) )U\(t>) (10) 

whcrc a« is defined by and the final step is true because XZ = 

luZX. Here we have succeeded in moving the generalised Pauli operators to 
where we want them; the only trouble now is that the overall gate that we 
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are implcmcnting is not identical to the one we want. We can get round this 
by using an adaptive measurement: Given x in advance, we can use one-dit 
teleportation to implcmcnt the gate FZ(aS~ x ') instead of FZ(a). With this 
adaptation, our new state is then given by (up to phasc) 

= X x 'z z 'FZ(a)U\<p) 

where 

x = m + z mod d; z — —x mod d. (11) 

So, by keeping track of the exponents of the X and Z operators by classical 
computation, and adapting the measurement taken via the original exponent 
vàlues, we can apply the gate FZ(&) to our state, up to a generalised Pauli error. 
The final step of the algorithm, which is (at least in all quantum algorithms 
to date) a measurement on the computational basis, can simply be performcd 
without any quantum corrections to the final state, because the Z errors have no 
effect on the probabilities of each outcome (they only change the computational 
basis vectors by a phase) and the X errors permute the computational basis 
elements, which can be compensated for by a final classical computation. 

The above process would also work if we had a multi-qudit state |^), with 
generalised Pauli errors acting on each qudit, and we entangled one of its qudits 
with a further | +) qudit via a controlled-Z interaction, and then performed one- 
dit teleportation. This would simply implement a one-qudit operation on the 
given qudit, and change the generalised Pauli errors on that qudit also. This 
situation often arises when using more complicated multi-dimensional clusters, 
and is illustrated by Figure|ü] 

We also note that we can move generalised Pauli errors through controlled-Z 
gates, because of the identities 

CZ(Z (g> I) = (Z ® T)CZ\ CZ(I ®Z) = (I® Z)CZ; 
CZ{X®I) = {X®Z- l )CZ; CZ{I®X) = {Z- Í ®X)CZ (12) 

which are easily established directly. These identities are necessary since we will 
sometimes need to treat a controlled-Z interaction as an actual two-qudit gate, 
rather than a means to allow the use of one-dit teleportation. An example of 
this (from Figure[5j) is the vertical interaction between rows (logical qudits): we 
use these to create two-qudit interacting gates. The above relations mean that if 
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Figure 6: One-dit teleportation applied to one qudit of a three-qudit state. It 
is easily seen (e.g. by writing I*) = Ey \i) |j) \ipij)) that one-dit teleportation 
applied to one qudit of the multi-qudit state gives us the appropriate transfor- 
mation on the transported qudit i.e. |*') = (/ <g> I ® X m FZ(a)) |\ff). 
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two logical qudits are in a state (X Xl Z Zl ®X X2 Z Z2 ) \^), applying a controlled-Z 
gate to this and commuting it past the error operators yields a state 

(X X1 Z Z1 - X3 ® X X2 Z Z2 - X1 )(CZ) |*) . (13) 
3.2 Adapative computation 

Given that a quantum algorithm ends in a mcasurcmcnt in the computational 
basis, generalised Pauli X or Z errors are not the only error operators that can 
be compensated for by a classical compuation. Any overall operator that only 
permutes computational basis elements (and possibly modifies them by a phase) 
is also correctable by a final classical computation. 

Let us develop this idea further. Suppose we wish to implement a singlo 
unitary gate U via a series of measurements on some cluster state, but end up 
with Z k P p U, where p e Sd and 

P p = J2\p(l)){l\. (14) 

1=0 

We can then compensate for this at the end of a quantum algorithm via classical 
computation also. 

With this in mind, we can ask ourselves if we can somehow adapt our scheme 
for computation to include correction for more general permutation errors of the 
form above. We can express this more formally as follows: suppose we are given 
a state \<j>) = Z n P p U \ip), and we wish to apply the gate FZ(a) to the state via 
one-dit teleportation. We wish to establish whether we can peform a one-dit 
teleportation to obtain the state \<p') = X m FZ(sL)Z n P p U \ip), and then rewrite 
this in the form \<j>') = Z m 'P p ,FZ(a)U \tp), with p' e S d . 

A simple class of permutations for which this is possible can be defined as 
follows. Let c is a unit in the ring Z^. 4 It then follows that the operator 

d-1 

s c = ]T|cfc>(fc| 

fc=0 

is a permutation operator. Furthcrmorc, since 

S c -iF = ^=Y J \c- 1 k){k\Y,^ ml \rn){l\ 



V l,m 
l.m 



A unit in a ring is an element that has a multiplicative inverse within the ring. is equal 
to the quotient Z/(dZ), and can be thought of as the ring of modulo d arithmctic. In Z^, the 
units are the integers coprime to d (i.e. have no common factor with d). Hence, when d is 
prime, all non-zero elements are units; for non-prime dimensions not all elements are units 
e.g. 2 in Z4. If c is a unit, then cj = ck if and only if j = k. This property establishes S c as 
a permutation operator. 
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it may be possible to implement the gate F c rather than F by using the identity 
F = S c S c -iF = S C F C , and treating the remaining S'c permutation as an error 
operator (in the same way we treat powers of X and Z). We use the phrase 
adaptive computation to refer to this idea of implementing a different quantm 
gate by changing the leading error permutation that we correct for classically. 

Let us suppose we are given a state \ip) = X X Z Z U \<j>), which is in a linear 
cluster, and we measure it to induce a one dit teleportation, so the state becomes 
íbv eauation lfTÜ|) ) — X m+z Z~ x FZ(aS x ^)U \<p) (if the measurement outcome 
is m). Note that the powers of the X and Z operators have already been 
updated. By applying the identity F = S C F C here, we obtain the state 

|^') = X m+z Z- x S c F c Z(a^)U \(/>) . (15) 

From this form we can see that by considering the S c operator as an error, our 
computational state is equal to F c Z(aS x ')U \(f>), modulo a permutation error 
operator of the form X X Z Z S C . To be able to maintain this form after later one- 
dit teleportations and further uses of the identity F = S c >F c i (c 1 any other unit 
in Zd), we need appropriate relations between S c and the operators F, Z(a) 
that appear through the measurement process; however, it can easily be verified 
that 

FS C = S e -iF; Z{a)S c = S c Z{a!) « = a ck ). (16) 

which allow us to move the S c operator to the left of the intended unitary gate 
(the change in a above will mean adaptive measurement is needed here), and 
furthermore since S c Sd = S ce i, it follows that we are always able to obtain an 
error operator of the form X X Z Z S C , where c = 1, . . . , d — 1 and a unit in Zg. 

Note that we are still using the same quantum process (one-dit teleporta- 
tion), but a different classical correction procedure. This means that we can 
implement the gates F^ Z(a) by choosing the measurement bases appropriately. 

We note that for qubits, since the only unit in Z2 is 1 = — 1, this means that 
S—i = I, and since F^ = F, this effect does not exist for cluster computation 
using qubits. 

With this success for single-qudit operations, however, comes a word of warn- 
ing for multi-qudit operations, through the identities 

CZ(S C <g> I) = (S c ® I)C[Z C ] ■ CZ{I® S c ) = {I® S C )C[Z C ] (17) 

and since in our model the interactions between qudits are fixed, these changes 
in the effective interaction between e.g. neighbouring linear clusters are an 
artefact of this framework. However, since these extra permutation operators 
are not introduced by a measurement process but instead through our choice 
in the classical computation, we can in theory design quantum algorithms to 
get round this problem or even utilise it to allow us to implement different 
interactions between neighbouring logical qudits. For example, since S_i = /, 
applying F' rather than F twice in a linear cluster produces two cancelling S-i 
factors which we no longer need worry about. 

The group of permutations generated by the shifts X 3 and the multiplication 
maps S c is of order < d(d — 1), which is in general less than dl, the order of 
the symmetric group Sd- This means that the above does not deal with the 
most general permutation error. One issue with trying to correct for general 
permutation errors is that we would have to find a commutation relation between 
P p and F for general p, which seems hard to ascertain. 
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We could have chosen to place the S c operator in front of the Pauli operators; 
in this case, we would need the following commutation relations: 

S C X = X C S C ; S C Z = Z C ^S C . (18) 

These relations mean that S'c (where c is a unit in Z<j) is a member of the 
generalised local Clifford group Cd, which we define by the normaliser of the 
generalised Pauli group Vd = {uj a X b Z c \ a, b, c € Z d } i.e. 

C d = {Ue U(d) | UPtf e V d V P e V d }. 

In this case, S c corresponds to U, and products of powers of X and Z correspond 
to P. By looking in this group, we may find more permutations that we can 
correct for. However, it can be shown that in prime dimensions the operators 
Z, X, F, S c and P defined by P : \j) — > u^Ü +1 )/ 2 \j) are sufhcient to generate Cd 
(see appendixQ . Since our permutations are already built out of products of X 
and S c , we cannot obtain any more single qudit permutations from the Clifford 
group. This strongly suggests (but does not prové) that we cannot correct for 
any of the other single qudit pcrmutation operators within cluster computation. 
It is harder to specify the Clifford group more generally, and so it may be the 
case that in non-prime dimensions one can implement further quantum gates 
by introducing other types of single qudit permutation error operators. 

3.2.1 Adaptive computation using multi-level permutations 

When working with multi-qudit quantum algorithms, we can similarly intro- 
duce multi-qudit permutation operators. Given an n qudit state EU \^}, where 
E represents an error operator on the n qudits, and U S U(d n ) the desired 
quantum evolution acting on an initial state l'í), we can write 

EU\$) = EP.P~ l U\^) 

i.e. we let P~ l become part of the quantum evolution, and P part of the error 
operator E. To allow us to maintain this form of error operator we need P to 
have appropriate relations with F and tensor products of Z(a). For example, 
consider a two qudit permutation P. After moving P left through (Z(&) <E> 
Z(h)), the operator to the right of P msut be in tensor product form, so the 
transformations can be implemented using single qudit measurements. Up to 
permutations on individual systems, the only non-trivial permutation that does 
this is the swap operator V, since V(A B) = (B ® A)V . This gives us a swap 
of two qudits for free at any time, and the above relation for V can be used to 
both update leading error operators and to calculate how measurement patterns 
need to be changed due to the introduction of the swap operator. However, we 
again have an issue with interacting gates between logical qudits changing form 
(for example, through the identity CZ^^Ytí = VisCZx,z] the indices represent 
the two systems each of these unitàries operates on). As before however, it 
could be that this is a help rather than a hinderance if utilised correctly. The 
swap can be introduced even in the qubit model, and, more generally, we can 
introduce any system swapping operator in this manner for an arbitrary number 
of qudits. 5 

5 However, at the end of a computation before the final computational basis measurement, 
any permutation can be introduced in the above manner and corrected for classically, since 
there is no need to propogate the intoduced operators through other operators. 
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3.3 Example - a variation on Deutsch-Josza 



Here we give an example of using one-dit teleportations to calculate the effect 
of a cluster state computation, and how adaptive computation can be useful in 
implementing quantum algorithms on qudit cluster states. 

Suppose we have a function / : l? d — > Z^; (x, y) i— ► (x — a)(y — 6), and 
we are interested in finding a and b. Classically, this will require at least two 
evaluations of /. In this section we will show that the quantum algorithm 
presented in Figure El can find (a, 6) with one use of the function /, and that it 
can be implemented on a cluster state of few qudits. 



\+d-i) - f 



F 1 



Z~ 



F 1 



\+d-i) 



Figure 7: The quantum circuit for finding a, b. 



Let Uf be the unitary gate corresponding to the controlled evaluation of / 
(the first gate in Figure EJ. It can easily be verified that 



U f \xx) \x 2 ) 



-d-lí 



,ab 



\ + d-l) 

-b\ 



= u ab {Z- a ®Z- b )CZ\x l )\x 2 )\+ d - X ) 



and using this result it is easy to verify that this quantum algorithm maps the 
first two qudits to |— a) \—b) (up to a phase), and so the final measurements 
yield the vàlues of (a, b). 

The form of the transformation induced by Uf on the first two qudits is 
ideal for implementation on a cluster state depicted in figure [SJ The two rows 
will form the upper two (logical) qudits, and the vertical interactions will form 
interacting gates. In place of the unitary evaluation of / will be a measurement 
on qudits 1 and 3, since we can write Z k = Z(z k ) with (z k )i — 2-ïïkl/d. We will 
hence measure these qudits in the basis defined by (FZ(z k )Y. To implement 
the F^ gates, we then measure in the F^ basis on qudits 2, 4, and use adaptive 
computation in the form F = S-iF-i = S-iF^ . Let be the result of the 
measurement on qudit k. We can then calculate the evolution of the cluster 
state algorithm (ignoring any overall phases), updating the error operators as 
we proceed using the relations established in the previous section: 



1 2 F 



3 4 F 

Figure 8: The cluster state used to implement our quantum algorithm, complete 
with numbering for qudits (see text). The qudits labelled F represent the final 
two qudits at the end of the computation. 



13 



cz 



mcasurc 1,3 



cz |+)|+) = i !+)!+*) 



-in 



mcasurc 2,4 



(X mi <8> X m3 )(FZ(z- a ) ® FZ(*-"))CZ |+) |+) 
(r 0+m ' ® X- 6+m ^_!)(F ® F+)-^= V |+) |+fe) 

(I-" +mi ® X- t * m *S- í )CZ |+) |+) 

czpsr a+mi ® x- b+m3 5_i)cz |+) |+) 

( X -a+rn lz b-rn 3 g X' 1 *™* Z a ~ mi S-l)C\Z- X ]C Z |+> |+) 

(X m2 F® X m4 F)(X- Q+mi Z fc ~ m3 ®X- fc+m3 Z a -" ll S'_ 1 ) |+) |- 

= |— 6 — m2 + mj, mod d) |— a — + mi mod d) 

The vàlues of a and ò can then be read off from the final qudits. We note 
two things about the above computation: 

1. The measurement of qudits 1 and 3 introduce extra F gates on both 
logical qudits not present in the original algorithm. However, because 
(F ® F^)CZ |+) |+) = CZ |+) |+), by using the identity F = S-iF* on 
one qudit, these extra gates do not affect the rest of the computation, 
except for the fact that the first qudit corresponds to the value of b rather 
than a as in the circuit of Figure (and vice versa) . 

2. The original algorithm requires the use of a controlled-Z _1 gate, where 
as our cluster state only has controlled-Z gate between the logical qudits. 
However, pulling through the S-\ error operator (introduced in the pre- 
vious step) to the left through the CZ gate between qudits 2 and 4, the 
CZ gate changes into a C\Z~ X \ gate. 

While this example thought of purely as a quantum algorithm is a bit of 
a toy example, it does illustrates not only how adaptive computation can be 
helpful in the implementation of algorithms on cluster states, but also how the 
non-trivial commutation rclation between S c and CZ can actually prové to be 
of hclp rather than a hinderance. 



4 The stabiliser formalism for cluster state QC 

Up to this point, we have presented a completely self-contained approach to 
cluster state computation. Most of the current papers ^ HJ EI] take a much 
more algebraic approach using the stabiliser formalism, developed by Gottesman 
|17j in the context of error-correcting quantum codes. The stabiliser formalism 
allows us to describe cluster states through sets of eigenvalue equations, and 
manipulation of these eigenvalue equations leads to a theorem that can be used 
to establish how a particular set of measurements can implement a quantum 
gate on a particular cluster state. In this section we will briefly summarise this 
approach, and compare the mèrits of the two. 
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4.1 Representation of cluster states using stabilisers 



Before we describe how stabilisers can be used to represent cluster states, let us 
give a more formal definition of cluster states. Let G — (V, E) be a graph; the 
set V are the vertices (tiris corresponds to a labelling of our physical qudits) 
and E Ç V x V are the edges of G (corresponding to which pairs of vertices 
have a controlled-Z interaction between them). For a 6 V, let |.) represent the 
state of qudit a. Then our cluster state 140 c(G) can ^ e wr itten as 



'C(G) 



n cz (a,b) 

, (a,b)eE 




where CZ< a ^\ represents the controlled-Z gate between qudits a and b. 

It can be shown Plll0| that \4>)c[G) IS un iiuely determined (up to phase) by 
the eigenvalue equations 



x\ <8 (g) z h \\i>) = m 

\beN(a) J 

where U a is the application of the unitary U to qudit a, and N(a) = {b G 
V | (a, b) S -E} represents the neighbours of qudit a i.e. the qudits b connected 

to a by an interaction (edge in G). The operators S(a) = X\ ® (®b e Ar( a ) ^fcj 

are the stabilisers of the state \4>) c(G)i anc ^ totaJly determine up to an irrelevant 
phase the cluster state. 



4.2 Measurement patterns in the stabiliser formalism 

Before we can state the main result that shows how stabilisers can be used to 
hclp find measurement patterns that implcmcnt useful quantum algorithms, we 
need to establish a few definitions. Given a cluster state \<t>)c{G)i ^ us break 
its graph up into three parts: an input cluster C/(G) and output cluster Co(G), 
and the remaining body of the cluster, which we will denote Cm(G). We will 
number the n qudits in C/(G) and Co{G) from 1 to n, where qudit i in the 
output cluster can be thought of as qudit i from the input cluster after the 
quantum computation has concluded. We suppose further that the n qudits in 
Ci{G) can, before any interactions between qudits is applied, be initialised in 
the more general entangled state |í/>(in)). 

We define a measurement pattern on a cluster C(G) to be a set of the form 
of a function : V — > U{d), such that the columns of the unitary matrix 

Mc(G){ a ) form the measurement basis on qudit a in the cluster. 6 If the mea- 
surement outcomes form a vector m, let P(Mc<g)i m ) represent the projector 
onto the states corresponding to the particular measurement outcomes. 

We are now ready to state the main theorem for cluster state QC and sta- 
bilisers: 



This definition difïcrs slighlty from that used in other papers but is completely equivalent. 
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Theorem 1 I10p. Suppose the state IV')c(G) = (^Cj(G) ® -P(-Mc(G)j m ) ® 
-^Co(G)) l0)c(G) °beys the 2n eigenvalue equations 

X l· ,C 1 (G)(UX l UÏ) Coi G)Wc{G) = ^IV')c(G) (19) 
Z !c í (G)( t/Zít/t )c·o(G) l^)c(G) = a;A *' 1 l^)c(G) ' 

(20) 

Then, if M-Ci(G) = {(h^i) I * £ CV(G)}, and í/ie measurement outcomes on 
Ci(G) are given by Si, then after all measurements on Cj(G) and Cm{G) of 

1$ C (G)> £/ie -Z" 1 ^ sía<e IV'(o«í)>cf ( G ) géuen 6y |-0(ouí)) Co(G) - C/Í7 e |V>("i)), 
w/iere 

n 

[/ e = 0(Z i ) A -- s *(X í ) A - i . 

i=l 

The unitary gate U and the numbers \ x .i, ^zà can all depend on the outcomes 
of the measurements on the input and body parts of the cluster state. Many 
examples of the usc of this formalism can be found in [21 EU • 

4.3 A comparison of the two approaches to cluster state 
QC 

The two approaches to cluster state QC that we have mentioned are very dif- 
ferent. The stabiliser formalism is a powerful tool in establishing cluster states 
and measurements on them to implemcnt particular quantum gates, but as ex- 
amples in 2 show, using it is an involved and difHcult process. The one-dit 
teleportation approach described here is a much more elementary description of 
the model, and it has two important advantages over the stabiliser approach. 
Since this approach is very closely linked to the circuit model for quantum com- 
putation, it means that it is possible to use some of the intuition and results 
from this line of thought in designing quantum algorithms using cluster states. 
This approach is also the easiest to understand as it uses no more than some 
bàsic linear àlgebra to establish the theory behind the model. 

The fact however that the one-dit teleportation approach is so closely linked 
to the ciruit model means that this approach is not suitable for considering 
all possible cluster states. Up to this point, we have always implicitly con- 
sidered clusters that are in a grid (two-dimensional), such that each row of 
cntanglcd qudits corresponds to a logical qudit, to which we apply successive 
one-dit teleportations, and any vertical interactions (columns) form two-qudit 
gates between logical qudits (as in Figure^J). However, clusters where a row 
does not form a logical qubit of the computation, but forms part of an inter- 
action between logical qudits, cannot have one-dit teleportation ideas applied 
to it. The cluster in Figure|Ü|is an example of this. The central qudit causes 
the problem: its presence means that there are not two disjoint paths from left 
to right through the cluster such that every qudit is in one path, so we cannot 
treat the cluster by implementing one-dit teleportations along the paths (to give 
single qudit operations) and considering links between the paths as interacting 
gates. 

Another issue with the one-dit teleportation approach is that of the timings 
of the measurements. In all of our clusters so far, all measurements in a row 
of entangled qudits must be taken 'from left to right' i.e. we must apply each 
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Figure 9: An example of a cluster that we cannot evaluate using one-dit tele- 
portations alone. 

teleportation in turn. While we could make joint measurements of qudits where 
no adaptive measurements take place, we cannot measure the qudits out of 
turn in the teleportation approach. Examples in [2] using the stabiliser method 
however distinctly allow measurements to be taken in an arbitrary order (up to 
the issue where some measurements may depend on other measurement results). 

Importantly though, the tool of adaptive computation does not naturally 
arise within the stabiliser formalism. It seems difficult to incorporate the extra 
operator errors that arise from adaptive computation into the major theorem 
above (particularly as their deterministic origin differs from the randomness of 
the origin of the Pauli operators), and it remains to be seen if there is an elegant 
or natural way that these two methods could be combined. The fact that the 
same quantum process is taking place but a different correction procedure is 
implemented would somehow have to be reflected in any such grander theorem. 

5 Universal gates 

In this section we are going to discuss how we can implement any unitary evo- 
lution on a single qudit using linear clusters of qudits. Universality for many 
qudits using multi-dimensional clusters follows from the comments in section 2. 

5.1 An approach to generating unitary gates 

In this section we are going to present a couple of bàsic results that present a 
possible method of generating unitary transformations on any QC system. 
We start by quoting the following result: 

Lemma 2 (p48, [19j). Let Hk be a basis over the reals for the space of d x d 
Hermitian matrices. Then any U G U(d) can be formed by products of the 
matrices exp(i/3kHk) for real /3k- 

Now, given a matrix U, with eigenvectors \Uj), we will define 



Suppose we can find a (minimal) set of matrices U% , . . . , U m such that the pro- 
jectors onto the eigenvectors of these matrices form (over the reals) a spanning 
set for the space of d x d hermitian matrices. It follows from lemma that 
products of the unitàries from the set 




(21) 



C = {U j (a)\j e{l,...,d}, ae[0,2if} 



17 



is enough to generate any unitary in U(d). 7 

This result in itself does not specifically apply to our model of QC yet; we 
need to add more constraints to the kinds of bases that we can use. 

5.2 Mutually unbiased bases and Hermitian matrices 

To motivate this discussion, we will state here (and prové later) that, using 
a linear cluster of three qudits and two measurements, with the first (left- 
most) qudit in a given state \ip), we can implement the gates Z(a) and X(a) = 
Y,k eia "\+k)(+k\ for any dimension d. The two bases {\k)}fzl and {\+k)}Í=o 
are what we call mutually unbiased, a property defined as such: 

Definition 1. Two bases {|«fc)}fcZj and {|frfc)}f=o are mutually unbiased if 

\(ai\bj)\ = 1/Vd for all 

Indeed, if we accept one-dit teleportation as the basis of the cluster state 
model, the only bases that we can measure are those that are mutually unbiased 
with respect to the Standard basis. We can take this connection further with 
the following observation: 

Lemma 3. Let {|^)}^ =1 be d+ 1 bases ofC d (with k = 1, . . . , d+1), such that 
they are mutually unbiased. Then real combinations of projectors of all these 
vectors span the space of d x d Hermitian matrices. 

This arises from the fact that mutually unbiased bases arise in the problcm 
of quantum state determination. In the limit of many measurements, measure- 
ments in each of a full set oï d+1 mutually unbiased bases are enough to specify 
any Hermitian matrix |20| . So mutually unbiased bases are potentially a route 
into finding a basis of projectors for the set of Hermitian matrices that are 
appropriate to our model. 

We will list briefly some known facts about mutually unbiased bases. It is 
known |2U1 I21| that the maximum number of possible mutally unbiased bases 
in dimension d is d+ 1. In dimensions that are prime powers, a construction is 
known that gives this maximal possible number of such bases [231 EU- However, 
more generally no construction is known to find such a maximal set in arbitrary 
dimension. For the first non-prime power case, d = 6, there are 3 known MUBs, 
but it is not known if any larger set exists. 

5.3 An explícit construction for prime dimensions 

For dimension d a prime, it can be shown via a result of |21| that a set of 
d+1 mutually unbiased bases can be obtained via the eigenvectors of the d+1 
matrices 

Z, X, zx, zx 2 ,... 1 zx d - 1 . 

The bases {\k)} d .ZQ and are the eigenvectors of Z and X respectively. 

What we will show in this section is that we can generate the gates Z(a), 
X(sl) and ZX k (a) for prime dimensions on linear elusters. Then, by the results 

7 Note that we not using operators of the form e.g. exp(iaZ) for some real a, as is used for 
qubits. This is because for d > 3, Z, X are not hermitian, and so exp(iaZ), exp(iaX) are not 
unitary matrices. 
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of the previous two sections, this will establish single-qudit univcrsality on linear 
clusters as required. 

It is easy to implement Z(a) and X(a) in any dimension using two measure- 
ments on a linear cluster of size three, because we may write 

Z(a) = F^.FZ(a); X(a) = FZ(a).F f . 

and hence by appropriate one-dit teleportations and use of adaptive compu- 
tation we can implement these gates. 8 We note that this scheme uses fewer 
measurements to implement these gates than the one given in These two 
gates will be the building blocks for implementing the gates ZX k (a). The fol- 
lowing lemma is crucial in allowing us to do this: 

Lemma 4. Let {\ip k )}j =1 be the eigenvectors of ZX k for k = 1, . . . ,p— 1, with 
eigenvalues lü ' j . There exists a phase transformation Z(b k ) such that 

Z(b k )\+ 3 ) = \^ k ). 

The index jk in the ket determines which eigenstatc of ZX k the phase 
transformation Z(b k ) maps \+j) onto; importantly, the index is not solely j, 
and depends on k. The proof of the lemma is a little technical and hence we 
postpone it to appendix[5] The result of this is that we may write 

ZX k ( a ) = Z(b k )X( a ^^)Z(b k y 

where a( X:fc ) is a permutation of the elements of a defined by a\ x = a k -n. 
Hence we can decompose ZX k (a) as 

ZX k (a) = F^.FZ(b k ).FZ(a^'^).F^Z(-bu) 

and so we can implement these gates using four measurements and use of adap- 
tive computation. 

5.4 The d = 3 case 

In this section we will give the above results in dimension 3 in more detail, to 
show explicitly which measurements are needed to implement each gate. 

In dimension 3, we have four primitive gates in our construction: Z(a), X(a), ZX(a) 
and ZX 2 (a). We need the eigenvectors of each of the matrices Z, X, ZX and 
ZX 2 . The computational basis elements are the eigenvectors of Z. The eigen- 
vectors of X are the columns of F i.e. 

l+o) = ^(|0) + |1> + |2)) 
l+i) - j= 3 (\0)+u;\l)+^\2)) 
l+2> = ^(|0)+^|1)+^|2)); 

the eigenvectors of ZX are 

\ZX ) = ^(|0) + |l)+c 2 |2)) 
\ZXt) = -L(^|0) + |l) + |2)) 
\ZX 2 ) = -^(|0)+c 2 |l) + |2)); 

8 Note that the gate can be implemented without the use of adaptive computation with 
three one-dit teleportations, since Ft = F 3 , but this takes up more resources. 
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and the eigenvectors of ZX 2 are 



\ZXl) = -^(|0)+ W |1> + |2)) 
\ZX\) = i (|0) + |l)+a;|2» 
\ZXl) = ^( W | ) + |l) + |2)). 



By utilising the proof of lemma^in appendixlBl íor indeed by direct verification) 
defining bi = (0,0,4tt/3) t and b 2 = (0, 2tt/3, 0) T , so that 



Z(bi) 

we have that 
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Z{h x )\+ ) = \ZX ) 
Z(b 1 )|+ 1 ) = w\ZXi) 
Z{h x ) |+ 2 ) = \ZX 2 ) 



and 



Z(b 2 )|+ > = \ZXl) 
Z(b 2 )|+ X ) = w 2 |ZX 2 2 ) 

Z(b 2 )|+ 2 > = |ZX?>. 

From these relations we can deduce that 

ZX{a) = Z(bi)X(a)Z(bi) t 



and 

where a' is dehned by 



ZX 2 {sl) = Z(b 2 )X(a')Z(b 2 ) t 



ai; o 2 = a 3 ; a 3 = a 2 



and using the measurement patterns in the previous subsection we can imple- 
ment these gates on linear clusters. 

The issue of adaptive computation in this dimension is simple, because only 
when c = —1 is S c unitary and not equal to the identity. So we can imple- 
ment the gates FZ(sl) and F* Z(a) in this dimension, and by doing so we only 
introduce one new kind of error, namely that introduced by S-i. 

Furthermore, it is easy to show that the shift permutation (312) and the 
transposition (23) generate the whole of the permutation group S3; these per- 
mutations correspond to X and S-i and so when d = 3 we can in fact correct 
for all possible permutation errors. 



5.5 Other cases 

As noted above, constructions for full sets of mutually unbiased bases exist in 
all prime power dimensions, but their constructions are not as simple as the one 
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listed above. In |21| . an explícit description for a set of mutually unbiased bases 
in d = 4 is given, and can be seen by inspection that a relationship similar to 
that for prime dimensions given by lemma J3J can be given. It remains to bc 
seen whether this relationship exists more generally. 

Furthcrmore, these mutually unbiased bases should be compatible with our 
scheme for generating universal gates in the cluster state model. Ideally one 
of the bases should coincide with the Fourier basis |Q, to fit in with the F 
transformation from one-dit teleportation. This is not the case in the example 
in |21| , which makes the task of generalising the ideas in this paper to non-prime 
dimensions more difhcult. 

It is clearly not the case however that non-existence of a full set of MUBs 
implies that our cluster state model is not universal. It is established in ^ü] that 
our given model is universal for all dimensions. They also prescribe a method 
which is equivalent to finding a set of projectors that is not only a basis for 
the set of Hermitian matrices, but also suitable for the cluster state model. It 
remains to be seen whether the model can be described more directly as in this 
paper. 

6 Conclusion 

In this paper we have developed a framework for quantum computation using 
cluster states with qudits as our bàsic physical resource, and while we have 
recovered all of the features of the qubit model, we have found that the extra 
degrees of freedom in using qudits allows us to control the quantum computa- 
tion using adaptive computation in a way that is not possible (for single qudit 
operations) when using qubits. The cases of prime dimension lend themselves 
particularly well to a simple description, and as a result it is hopeful that these 
ideas could be taken further and used to design quantum algorithms on qudits. 

One important issue we have not covered in this paper is that of fault toler- 
ance. There are a number of works on this issue specifically pertaining to the 
qubit cluster state model [El ESI EH] > and given the similarities between this and 
the qudit model, it seems likely that a number of the results may well jump 
across to the qudit model. How the issue of adaptive computation will fit into 
the picture could however be a sticking point for generalising these results. 

The most technologically advanced physical implementation of cluster state 
QC is an optical solution as proposed by Nielsen 0] and refined by Browne and 
Rudolph [22] ■ A photon can be thought of as a qubit by treating the polarisation 
of the photon as the relevant degrees of freedom, and in |23] Knill, Laflamme 
and Milburn presented a scheme for entangling photons using measurement. 
This idea is at the center of Nielsen's proposal; Browne and Rudolph present 
another optical scheme that is more efficient in its use of resources (linear opti- 
cal elements). 9 However, by using other degrees of freedom present for photons, 
they could be potentially thought of as qudits. It would most definitely be in- 
teresting if one could find a scheme that is similar to any of the above proposals, 
but within which photons are treated as qudits, and furthermore, whether this 
idea lends itself to a scheme for qudit cluster state quantum computation. 

A number of other interesting qüestions arise from this model. One immedi- 
ate question is whether the computational paradigm that arises from adaptive 

9 Experimental work using photonic cluster states has been reported in 0. 
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computation lcnds itsclf to solving a particular class of problems. Furthermorc, 
we have not explored the potential for designing clusters with controlled-Z fc 
interactions, as was mentioned in section 3. Finally, we have not made any 
attempt here to consider fundamentally the role of entanglcmcnt within the 
cluster state model. This encompasses a large number of potential qüestions, 
and further investigation could lead to some insight into the relationship be- 
tween computation and entanglement. 
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A Maximal connectedness of cluster states 

This appendix establishes the result that if any one qudit in a cluster is mea- 
sured, the remaining qudits form a cluster state up to some Pauli Z errors. 
As mentioned in the body of the text, for cluster of qubits this issue was first 
discussed in ^^1; nere we gi ye a different but elementary proof of the same fact. 

Let \tp) be a cluster state. Since all the controlled-Z interactions between 
pairs of qudits commute, we can write this as 

\Ï>)=S\+)\tP') 

where is the cluster state of all but one of the qudits, and S represents the 
controlled-Z interactions with the remaining |+) qudit. We can write the above 
as 

d— 1 

= —=^\k)(z k ®...®z k ®i®...®i)W) 

where the Z k operators act only on the qudits that are immediately connected 
to the final qudit (that is, those that have been made to interact with the final 
qudit by a controlled-Z operation). It is clear from this form that measuring 
this qudit in the computational basis and obtaining result j will project the rest 
of the state onto a cluster state with Pauli errors on all qudits that were 
adjoined to the measured qudit in the original cluster. 

Since the cluster state of two qudits is equivalent (up to unitàries on the in- 
dividual qudits) to a Bell state, our notion of maximal connectedness is identical 
to that discussed for linear clusters in |12j . 
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B Proof of lemma H] 

To prové the lemma, we first need to compute the eigenvectors of ZX k in prime 
dimensions d > 3. 

Lemma 5. Define \a^) by 

d-i 



1=0 



where a satisfi.es ai+k + 1 = oli mod d (with mod d arühmetic in indices also). 
Then |a«) is an eigenvector of ZX k with eigenvalue 1. 



Proof. Rewrite |a^) as 



(k)\ = J_ 



QL / 



■iy/'+' \i+k) 



(with modulo d addition as appropriate). Then 



d-i 



ZX k 



d-1 



The other eigenvectors of ZX k are then given by X m \ç/- k ^), since 



zx k .x- r 



W> = uj m X- m .ZX h 



m xa — m 



v(*)\ 



For definiteness, we define the phases of the eigenvectors by letting \4>o) = 

|a^), with a k ^ — 0, and \ip k ) — X~i \g^ k ^)- The following lemma is equivalent 
to lemma 0] 

Lemma 6. Let {\ip k )}^Zo be the eigenvectors of ZX k for k = 1, . . . ,p — 1 as 
defined above. Then there exists a phase transformation Z(bfc) such thai 



Z(h k ) 



LO 2 



Proof. By the above lemma, 

d-l 



V " 1=0 
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where a>i + k +1 = ai mod d. Define by (bk)i — ai. Then 

However, it is easily seen (e.g. by induction) that, for any r, s, 

at = ai +jk +jl+~j(j-l)k (22) 



and hence 



i=0 

1 d-l 

d-1 



v/d _ 

1 



í=0 



So we have that 



z(h k )\+ 3 ) = u J y^ k \^ k ) 



and since the dimension here is prime, as j runs over all vàlues 0, . . . , k — 1, 
so does j/c. We can recover lemma 0] from this by the phase transformation 

|^*) i— > ui^ J ^ k |V$}> hence removing the phase that appears in the above 
lemma. 

The eigenvectors of ZX for d = 2 are given by |-0±) = ^ (|0) ±i|l)). 
The above proof fails because (|22J) is only self-consistent when d is odd: By 
substituting j for d, then modulo d, we obtain that ai = ai+dk+dl+^d(d—l)k = 
ai + nd(d — í)k mod d i.e. we require ^d(d — l)k = mod <i, which is true 
if and only if d is odd. However, lemma 0] still holds: let b = (0, tt). Then 
Z(b) = |0)(0| + ï|l)(l|, and it is easy to verify that 

Z(b) |+o) = |^+> ; Z(b) = \ï>-) . 

C The Clifford group in prime dimension 

The sole result of this appendix is a result of Gottesman JHj : 

Lemma 7. The Clifford group C'd (normaliser of the generalised Pauli group 
T'd) is, when d is a prime, generated by the operators Z, X, F, S c and P defined 
byP:\j)^^+D/*\j). ' 
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Proof. The one-qudit Pauli group consists of elements of the form ui a X h Z c . Let 
us (for now) ignore the preceding powers of lo. Then, under conjugation by a 
Clifford group operation (i.e. for U € Cd, the transformation is the conjugation 
P i ► UPW), wc must have that (up to phascs) X ^ X l Z^ and Z ^ X k Z l . 
For A, B e Vd, dcfinc a(A, B) by AB = cj a(A ' B ^BA. Undcr conjugation by a 
Clifford operation, this commutation relation is preserved and so we must have 
that a(X i Zi,X k Z l ) = a(X,Z) = 1. Furthcmore, 

a{X i Z j ,X k Z l ) = a{Z\X k )+a(X\Z l ) 

= jka(Z, X) + ila(X, Z) = il - jk 

and so a general Clifford operation is an element of Vd (to choose the phases 
for the images of X and Z) times some operation from some class of operations 
labelled by k, l), with il — jk = 1. We know that the Pauli operators are 
in the Clifford group, and so we can ignore any phases as the Pauli operators 
can correct for these. We split up considering the full set of operations into 
two groups. First let us assume that i =/= 0. ft can be shown via elementary 
methods that if we define C(i,m,n) = S t P m Q n (where Q = FPF^), then, 
under conjugation by C(i,m,n), 

XHl*r lm ; Z i-> x in Z i ~ 1{1 - mn "> 

(note that sincc i ^ 0, i always has an inverse). By choosing m = — ij and 
n = í- l k, wc find that X ^ X' L Z^ and Z ^ X k Z\ with l = i' 1 ^ + jk) 
as rcquired. This gives us all Clifford operations for i ^ 0. For i = 0, we 
must then have jk — 1, and so j ^ 0. We wish to implement the operation 
X i ► Zi and Z X k Z l . We can do this by first performing X X~i 
and Z i ► Z k X~ l (which can be done from the calculations above) and then 
conjugating by F to implement the required transformation. This gives us all 
Cifford group operations. 

The above proof does not work for non-prime dimensions because of the need 
for inverses of all non-zero elements in Z^, which only exist if d is prime. 
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